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Abstract 

The fluid mechanics of multiphase flow breakup of 
liquid nitrogen, LN 2 , jets injected into sonic velocity 
nitrogen gasflow, was experimentally investigated. A 
scattered-light scanning instrument was used to measure 
the characteristic dropsize, D y 5 , of LN 2 sprays and to 
determine the effect of droplet vaporization on experi- 
mental dropsize measurements. Under sonic gas-velocity 
conditions, liquid-jet breakup occurred in the regime of 
aerodynamic stripping. As a result, the following 
correlation of volume-median drop diameter, D y 5 , with 
atomizing gas flowrate, W g , was derived for two-fluid 

atomizers; D y * = k c W g , where proportionality 

constant k c and exponent n are functions of droplet 
vaporization rate. In this study, partially vaporized 
sprays were investigated and it was found that 
n = 1.11, which is considerably less than the value of 
1.33 that is predicted by atomization theory. This was 
attributed to the evaporative loss of very small droplets. 
As a result, the following expression was obtained exper- 
imentally: D y * = 301 W* 11 . Values of D y 5 , that 

existed prior to partial vaporization of the LN 2 sprays, 
were calculated and the following expression was derived 

— 1 1.33 

for originally unvaporized LN 2 sprays: D y 5 = 285 W ’ . 

This expression agrees well with atomization theory that 
predicts n = 1.33, for liquid jet breakup in high-velocity 
gasflow. 

Nomenclature 

a acceleration, cm/sec 2 

drag coefficient 

D y 5 volume median drop diameter 
k c correlation coefficient, sec/g-cm 

m mass of characteristic LN 2 droplet, g 

n exponent for weight flow of fluid 
Nu heat transfer Nusselt number, based on D y 5 
Re Reynolds number based on D y 5 
t vaporization time, sec 

V fluid velocity, cm/sec 
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W weight flow of fluid, g/sec 

x axial downstream spray sampling distance 

7 kinematic viscosity, cm 2 /sec 

p fluid density, g/cm 3 

Subscripts 

a acoustic 

c calculated 

e experimental 

g gaseous nitrogen, GN 2 

1 liquid nitrogen, LN 2 

Introduction 

Two-phase flows of liquid nitrogen, LN 2 , and 
gaseous nitrogen, GN 2 , were experimentally investi- 
gated. LN 2 jets were injected into high-velocity GN 2 
streams inside a two-fluid fuel nozzle. The jets rapidly 
disintegrated into clouds of very small droplets. As a 
result, cryogenic sprays having relatively large droplet 
surface-areas per unit volume of liquid were produced. 
High surface- area fuel sprays are desirable since they 
provide rapid vaporization and efficient combustion of 
liquid fuel sprays in gas-turbine and rocket combustors. 
However, in order to calculate vaporization and burning 
rates of fuel sprays, dropsize data are needed that will 
accurately represent sprays in terms of a characteristic 
dropsize, such as the volume-median drop diameter, 
D y 5 , used in this study. Dropsize data were obtained 
for cryogenic sprays and used to derive correlating 
expressions for two-fluid atomization, in terms of 
interacting aerodynamic and liquid surface forces. Such 
information is needed for analytical modeling of fuel- 
spray combustion processes. Also, fuel-spray dropsize 
correlations are needed that will cover wide ranges of 
properties of liquid fuels and atomizing gases. This will 
aid in the development of more efficient fuel nozzles and 
thereby improve combustor performance. 

Highly volatile cryogenic liquid sprays containing 
relatively large numbers of very small droplets vaporize 
quite rapidly. This is due to the presence of high 
thermal gradients that exist in the gas-film surrounding 
the droplets. In the present study, the atomizing-gas 
was close to room temperature, 293 K. As a result, the 



original characteristic dropsize initially formed at the 
atomizer orifice was not measured directly. However, it 
was possible to determine the original value of D y 5 , by 
calculating the amount it changed due to evaporation. 
Vaporization rates were calculated from heat-transfer 
and drag coefficient data that had been previously 
obtained for drops accelerating and vaporizing in high- 
velocity gasflows. Such data are difficult to obtain, 
especially when the droplets are microscopic in size and 
attain relatively high velocity in a very short distance. 

In the present study of cryogenic sprays, character- 
istic dropsize measurements were obtained with a 
scattered-light scanner. It was originally developed at 
NASA Lewis Research Center. The instrument was 
improved optically, as described in Ref. 2. As a result, 
accurate dropsize measurements were obtained even in 
the presence of high velocity and high gas-density 
gradients. Therefore it was possible to determine the 
effect of droplet vaporization on experimental dropsize 
measurements, with the improved instrument. Also, it 
should be noted that sonic gas-velocity conditions were 
used in this study and liquid-jet breakup occurred in the 
regime of aerodynamic stripping. 

The effect of atomizing gas flowrate on liquid jet 
breakup in gas streams has been reported in Refs. 1 to 7. 
In Ref. 1, water sprays were formed with a two fluid 
atomizer. The volume median drop diameter, D y was 
measured for each spray and correlated with atomizing 
gas flowrate to give the following general expression: 

D ! = k r W n , where W„ is GN 0 flowrate. Atomiza- 

tion theory, Ref. 8, predicts that n = 1.83 for the case 
of liquid jet breakup in high velocity gasflow, i.e., in the 
regime of aerodynamic stripping. 

Numerous investigators, as indicated in Refs. 1 to 7, 
have studied the effect of gas velocity, V , on liquid-jet 
breakup in gas streams. Their results are summarized 
in Table 1, In the study of water-jet breakup that is 
described in Ref. 1, good agreement was obtained with 
atomization theory when the exponent n was deter- 
mined. However, in the case of LN 2 sprays, good 
agreement with atomization theory may not occur. This 
can be attributed to the fact that highly volatile sprays 
can become partially vaporized. In Ref. 1, the study of 
water sprays showed a marked effect of droplet vaporiza- 
tion on measurements of D y 5 , when sampling distance 
downstream of the atomizer orifice was increased from 
2.2 to 4.4 and 6.7 cm. 

In the present study, original values of D y 5 for 
LN 2 sprays, i.e., prior to their partial vaporization, were 
calculated over a range of atomizing- gas flowrates. To 
do this, heat transfer and drag coefficient expressions 


obtained from Ref. 9 were used to calculate changes in 
dropsize that occurred due to partial vaporization of the 
sprays. Thus, values of AD y 5 produced by droplet 
vaporization could be calculated and used to determine 
the original value of D y 5 that existed before vaporiza- 
tion occurred. LN 2 sprays were injected into a low 
velocity, 5 m/sec airflow. Fluid pressures for LN 2 and 
GN 2 were varied over a range of 0.2 to 1.0 MPa and the 
entire spray was sampled with a 4.4- by 1.9-cm rectan- 
gular laser beam. The center line of the laser beam was 
located at an axial distance of 1.2 cm downstream of the 
fuel nozzle orifice. 

Apparatus and Procedure 

As shown in Fig. 1, a two-fluid fuel nozzle was 
mounted at the center line of the 24-cm diameter duct 
and operated over pressure ranges of 0.2 to 1.0 MPa for 
both liquid nitrogen, LN 2 , and the atomizing gas, 
gaseous nitrogen, GN 2 . LN 2 sprays were injected 
downstream into the airflow. The two-fluid nozzle was 
fabricated as shown in Fig. 2, with an orifice diameter 
of 0.66 cm. LN 2 at a temperature of 77 K was axially 
injected into the airstream by gradually opening the 
control valve until the desired flowrate of 51 g/s was 
obtained as indicated by a turbine flowmeter. GN 2 was 
then turned on and weight flowrate was measured with 
a 0.51-cm diameter sharp-edge orifice. After the air; 
atomizing-gas and LN 2 flowrates were set, the volume 
median diameter of the spray, D y 5 , was measured at a 
distance of 1.5 cm downstream of the atomizer orifice. 

The scattered-light scanner consisted of a laser 
beam expander with spatial filter, rotating scanning slit 
and detector, as shown in Fig. 1. The instrument meas- 
ured light as a function of light scatter angle, by repeat- 
edly sweeping a variable-length slit in the focal plane of 
the collecting lens. The data obtained is scattered-light 
energy, which is a function of the scatter angle relative 
to the laser beam axis. This particle size measuring 
method is similar to that described in Ref. 10. 

Measurements of scattered-light energy normalized 
by the maximum energy and plotted against scattering 
angle can be used to determine the volume median 
diameter, D y 5 , as discussed in Ref. 11 and illustrated in 
Fig. 3. Dropsize distribution dispersion can also be 
determined from this plot. Also, according to this 
method, measurements are independent of particle distri- 
bution function. For a typical measurement, the scan is 
repeated 60 times per second to average out temporal 
variations in the energy curve. Thus, spray pattern 
effects were minimized by measuring an average charac- 
teristic drop diameter, D y 5 , for the entire cloud of 
droplets. 
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Calibration was accomplished with five sets of 
monosized polystyrene spheres having diameters of 8, 12, 
25, 50, and 100 /xm. LN 2 sprays were sampled fairly 
close to the atomizer orifice and they contained relative- 
ly high number-density droplets. As a result, it was 
necessary to correct for multiple scattering as described 
in Ref. 12. Measurement corrections also included Mie 
scattering theory when very small droplet diameters, i.e., 
<10 /im, were measured. Dropsize measurements agreed 
within ±5 percent, as indicated by reproducibility tests. 
Corrections for background scattering caused by high 
velocity gradients in the laser beam were made by 
obtaining new background readings, whenever gas 
velocity was changed and prior to turning on the spray. 

Experimental Results 

LN 2 spray measurements of experimental dropsize 
D y 5e were obtained with the scattered-light scanner 
located 1.5 cm axially downstream of the two- fluid fuel 
nozzle. Drop velocity, V d , was based on D y 5e and 
drag coefficients were calculated from expressions 
derived in previous studies using a high-speed droplet 
tracking camera developed at NACA Lewis laboratory. 
The technique is described in Ref. 9. Resulting values 
of V d and C d were used to calculate vaporization 
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time At, which was then used to determine AD .. 

v.o 

Finally, calculated values of D y 5c were obtained and 
used to write an expression for the originally 
unvaporized LN 2 spray in terms of the effect of atomiz- 
ing gas flow on D y 5c . 

Correlation of D y 5c with GN 2 Flowrate 

The entire LN 2 spray cross section was sampled at 
1.5 cm axially downstream of the two fluid fuel nozzle. 

Values of the reciprocal volume median diameter, D y 

were measured and plotted against GN 2 flowrate, W , 
as shown in Fig. 4. Since the tests were conducted in 
high velocity gasflow, LN 2 jet breakup occurred primar- 
ily in the regime of aerodynamic stripping. From the 
plot shown in Fig. 5, the following correlation of recipro- 
cal Dy.Se atomizing gas flowrate, W g , was 

obtained: 

= 301 W 111 ' (!) 

v.5e g 

where D y 5e and W g are expressed as cm -1 and 
g/sec, respectively. 

The exponent 1.11 for W g is considerably less than 
1.33, which is predicted by theory for liquid jet breakup 
in high velocity gasflow. This discrepancy is assumed to 
be due to a loss of small vaporizing LN 2 droplets, prior 
to spray measurements. In Ref. 2 it was found that at 


a LN 2 flowrate of 76.5 g/sec the exponent n agreed 
with atomization theory. However in the present study, 
the LN 2 flowrate was set at a lower value, i.e., 51 g/sec. 
This was due to the fact that high LN 2 flowrate readings 
often indicated excessive N 2 vapor entrainment in the 
LN 2 jets, which is detrimental to the atomization 
process because a high pressure drop and flowrate can 
cause “flash” boiling breakup of the liquid jet to occur 
at the nozzle orifice. This phenomena was not studied. 

Experimental results obtained in the present study 
agree better with atomization theory than those reported 
in Ref. 2. This is due to the allowance made in the 
present study for the effect of droplet vaporization on 
volatile liquid dropsize measurements. This effect was 
not accounted for in Ref. 2. As a result, dropsize data 
did not give agreement with the theoretically predicted 

gasflow exponent, n, in the expression, D y * = k W”. 

However, the proportionality constant k was too low 
to adequately characterize the original spray. This was 
attributed to the fact that the study in Ref. 2 did not 
take into account the effect of small droplets vaporiz- 
ing completely, before dropsize measurements could be 
taken. 

LN p Droplet Velocity Calculations 

To determine the effect of droplet vaporization rate 
on experimental values of D y 5e , it was necessary to 
calculate vaporization time, t, by First calculating 
droplet velocity, V d , for D y 5e over a distance of 
0.13 cm in the nozzle throat and a distance of 2.2 cm 
from the nozzle orifice, x = 0, to the downstream edge 
of the laser beam, as shown in Fig. 1. In order to deter- 
mine volume-median drop velocity, V d , the acceleration, 
a, of LN 2 droplets was calculated from the following 
momentum balance: 

mj. = i,.A d (v, - v 4 !c d (2) 

where m d and A d are mass and area of dropsize 

D v.5e> res P ectivel y> ie > m d = ^l ffD v.5t/ 6 ’ and C d is 
drag coefficient for D y 5e . Rewriting this expression, in 

2 

terms of change in drop velocity squared, A V d , over the 
distance of travel, Ax, the following relationship is 
obtained; 



3 



where C d = 27Re 0 * 84 , as given in Ref. 9, and Re is 
based on the characteristic dropsize, D y g . 

To determine momentum transferred from GN 2 flow 
to the LN 2 droplets, the deceleration of an atomizing 
gaseous nitrogen jet into a surrounding low velocity 
airflow was determined as follows. At the nozzle orifice, 
GN 2 velocity, V g was equal to the acoustic velocity, V c , 
of gaseous nitrogen. Values of V g that were used to 
solve Eq. (S), were calculated at downstream distances 
of x = 0, 2.2, and 5 cm, respectively and plotted in 

2 

Fig. 5. The calculated values of are based on data 

given in Ref, IS. These data are also plotted in Fig. 5 
for comparison. The percent deceleration of the atomiz- 
ing nitrogen gas is assumed to be approximately the 
same in both cases, since the two-fluid nozzles used in 
Ref, 13 and the present study were very similar in 
design. 

2 

By means of graphical integration, values of V d 

were calculated and plotted against downstream dis- 
tance, x, as shown in Fig. 6. LN 2 droplet vaporization 
time At was calculated from this plot and the expres- 
sion At = x/V d . Calculated values of At for a given 
distance Ax are given in Table 2 along with calculated 
Reynolds numbers and drag coefficients averaged over 
the Ax distance. Gas and liquid transport properties 
used in calculating vaporization times are given in 
Table 3. Also, it should be noted that initially, in the 
nozzle throat, i.e., Ax = —0.127 to 0 cm, the reciprocal 
value of D y g is the same as that given in Fig. 4 for 
Dy 5c . Whereas at Ax = 0 to 2.2 cm, it decreased to 
the value of D y ge . 

Spray Vaporization Rate Calculations 

The vaporization rate of cryogenic sprays charac- 
terized by D y g were obtained by using the following 
heat-balance equation: dm d /dt = hA AT/H t , where h 
is the heat-transfer coefficient, A is droplet surface-area, 
based on D y g , AT = T g ~ T d , and H fc — Hv 
■f Cp AT, where Hy is the latent heat of vaporization 
of liquid nitrogen and C p is the specific heat of gaseous 
nitrogen. This expression may be rewritten as follows to 
obtain vaporization rate in terms of changes in droplet 
surface-area with time: 

dD* 5 / dt = 4k g AT Nu / Pl H t ( 4 ) 

where k g and p j are gas thermal conductivity and 
liquid density, respectively. In a previous fuel-droplet 
study using a high-speed droplet tracking camera to 
determine vaporization rates for fuels such as n- octane 
and methanol, as well as numerous other liquids includ- 


ing water, benzene, acetone, nitrobenzene, butanol, and 
carbon tetrachloride, it was found that: Nu = 2 

4- 0.303 Re 0 ' 6 where Re = p g D v g AV/p^, and AV 
is the average velocity difference over an incremental 
distance Ax. 

To calculate vaporization time, t, inside the two- 

2 

fluid nozzle, an average value of V d was taken 

between x = —0.13 and 0 cm, and used to calculate 
At, i.e., At = Ax/V*. In a similar manner, the average 
Reynolds number, ne, was determined by averaging 
Reynolds numbers over a distance of 0.13 cm. Re and 
t were then used in heat-transfer calculations to deter- 
mine LN 2 droplet vaporization rates. Calculations were 
based on characteristic diameters D y gc , D y ge and 
transport properties listed in Table 3. GN 2 viscosity 
and thermal conductivity were evaluated at the average 
gas-film temperature, i.e., Tj = l/2(T g + Tj). Droplet 
surface temperatures were assumed to be near the boil- 
ing point, 77 K, while droplet sprays were being accel- 
erated and partially vaporized. The latent heat of 
vaporization of LN 2 was evaluated at 77 K and the spe- 
cific heat of nitrogen vapor was evaluated at the average 
gas-film temperature, T f . 

Calculated values of the change in characteristic 

droplet diameter squared, D y g , that occurs due to par- 
tial vaporization of the sprays, are given in Table 4. 
The calculation was made for gaseous- nitrogen and 
liquid- nitrogen flowrates of 1.82 and 51 g/sec, respec- 

tively. Calculated values of D y g , were obtained from 
2 2 2 

the expression; — AD yg = D y gc — D yS ^. Values of 

D y ge are plotted against GN 2 flowrates as shown in 

Fig. 4. From this plot, the following correlating expres- 
sion is obtained for calculated values of D y gc ; 

D~* = 285 W 1 - 33 ( 5 ) 

v.5c g 

In terms of the Sauter mean diameter, D 32 , it was found 

“1 1 33 

that Eq. (5) can be rewritten as; D 32c = 422 W ' 

Comparison of Eq. (5) with Eq. (1), shows that the 
exponent 1.33 for gas flowrate in Eq. (5) is considerably 
greater than 1.11, as given in Eq. (l) for experimental 
values of D y ge . Equation (5) is in good agreement 
with atomization theory and the proportionality con- 
stant, 285 is only slightly less than the value of 301 
given in Eq. (1). Thus, droplet vaporization had con- 
siderable effect on the exponent, n, but only a small 
effect on the constant, k. Also, it should be noted that 
Eq. (5) is useful in establishing the initial values of D y g 
that are needed in analytical modeling of cryogenic 
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liquid spray vaporization and combustion, on a macro- 
scopic scale. 

The expression, D y * = k c used in the present 

study of LN 2 sprays cannot be directly compared with 
expressions obtained in previous studies of water sprays, 
until the effect of atomizing-gas temperature on LN 2 
spray formation has been investigated and determined. 
In the case of LN 2 and water sprays, the exponent n 
should equal 1.33 in agreement with atomization theory. 
However, the value of k will be considerably different 
due to large differences in their properties. The effect of 
atomizing-gas properties on LN 2 jet breakup also needs 
to be further investigated. 

Concluding Remarks 

It was possible to determine the effect of droplet 
vaporization on cryogenic spray measurements because 
previous studies in Ref. 9 were able to provide heat and 
momentum transfer expressions for droplets accelerating 
and vaporizing in high velocity gasflows. As a result, an 
originally unvaporized volume-mean diameter, D y 5c , 
could be calculated for LN 2 sprays produced with two- 
fluid fuel nozzles. The results agree well with atomiza- 
tion theory and can be readily applied to analytical 
modeling of fuel spray combustion in gas turbine and 
rocket combustors, under the conditions used in the 
present study. Dropsize data for different atomizing 
gases and wide ranges of gas temperature are needed in 
order to generalize the expressions derived in this study 
of GN 2 and LN 2 . 
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TABLE 1.— ATOMIZING-GAS VELOC- 


TABLE 3.— ATOMIZING-GAS AND 


ITY EXPONENT, n, FOR LIQUID- 
NITROGEN JET BREAKUP, 


Dq ~ V 


Source (Ref.) 

Exponent 

Kim and Marshall (3) 

1.33 

Lorenzetto and Lefebvre (4) 

1.00 

Nukiyama and Tanasawa (5) 

1.14 

Weiss and Worsham (0) 

1.33 

Wolf and Anderson (7) 

a 1.33 

Adelberg, (theory) (8) 

1.33 

Present study 

1.33 


“Wax sphere dropsize data. 


TABLE 2.— CALCULATED LN 2 DROPLET 
VAPORIZATION TIME, AT, FOR D v 5 , 
AT W = 1.82 g/sec 

, ... 


Zone 

Ax, 

cm 

Re* 

C d 

AtxlO 4 , 

sec 

Inside nozzle 

0.13 

70 

0.70 

0.47 

throat 





Downstream 

2.2 

35 

1.36 

3.26 

of nozzle 





Total 

2.33 

— 

— 

3.73 


a Ke = RE X Afc 1 /t t + Re 2 At 2 /t t . 
^Inside fuel nozzle. 
c Outside of fuel nozzle. 


LIQUID-JET TRANSPORT PROP- 
ERTIES, AT W g = 4,54 g/sec, 

T g = 293 K, Tj = 77 K, 

AND x = —0.127 cm 

Nitrogen gas 
V , cm/sec . 

Pg. g/cc • • • 
a #*g> g/ cm sec 
8 k g , cal/sec si 

LN 2 droplets 
V^, cm/sec . 
p v g/cc . . . 

H,, cal/g .. 

C„y cal/g °C 

“Evaluated at T, = l/2(T g — T,). 

TABLE 4.— REDUCTION IN DROPSIZE 
DUE TO VAPORIZATION, -AD* 5 , 
AND CALCULATED ORIGINAL 

DROPSIZE, D^ “ 


g/sec 

-AD* xlO 8 , 
cm z 

D 2 v 10 8 
cnv 

Jc> 

cm 

1.82 

37.0 

289 

630 

2.72 

21.2 

113 

1044 

3.63 

19.3 

62 

1530 

4.54 

17.6 

39 

2162 


an” 1 — m2 _ AT)2 \ — 0.5 

u v.5c - ^v.Se 


3.43 xlO 4 

3.84 xlO 4 

1.25 xlO 4 

cm (°C/cm) 4.20xl0 5 


204 

0.80 

47.8 

0.25 
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Figure 1 . — Atmospheric pressure test section and optical path of scattered- tight scanner. 


Inside diameter 



0.663 cm 


Figure 2. — Diagram of pneumatic two-fluid atomizer. 


0 .005 .010 .015 

Light scatter angle, rad 

Figure 3.— Typical scattered -light energy curve 
with four characteristic diameters. 
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